STM imaging of impurity resonances on Bi2Se3 
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In this paper we present detailed study of the density of states near defects in Bi 2 Se3. In particular, 
we present data on the commonly found triangular defects in this system. While we do not find 
any measurable quasiparticle scattering interference effects, we do find localized resonances, which 
can be well fitted by theory 1] once the potential is taken to be extended to properly account for 
the observed defects. The data together with the fits confirm that while the local density of states 
around the Dirac point of the electronic spectrum at the surface is significantly disrupted near the 
impurity by the creation of low-energy resonance state, the Dirac point is not locally destroyed. We 
discuss our results in terms of the expected protected surface state of topological insulators. 

PACS numbers: 71.18.+y, 71.20.Nr, 79.60.-i 



Probing the interaction of quasiparticles with impu- 
rities and defects has proven to be a powerful tool for 
the study of electronic structure in solids. First, quasi- 
particle scattering from defects and impurities give rise 
to interference effects of electron-waves associated with 
bands dispersion of the solid. In addition, impurities and 
defects can also result in bound states that induce lo- 
calized resonances. For surface state bands (SSB), both 
phenomena are most suitably studied using scanningtun- 
neling microscopy (STM) and spectroscopy (STS) [2|-[4j. 
In the case of topological insulators (TI) [5, 6], probing 
quasiparticle scattering interference (QPI) and resonance 
scattering phenomena are expected to be particularly re- 
vealing, given that our focus is on the unusual properties 
expected from the SSB of these materials. 

Much of the recent STM work on impurity and defect 
scattering in TI systems concentrated on Bii_ x Sb^ Q 
and Bi2Te3 [1, [§]. Indeed, within the main part of the 
SSB, exhibiting convex contour, oscillations were shown 
to be strongly damped, supporting the hypothesis of 
topological protection. However, oscillations do appear 
in the upper part of the SSB exhibiting concave contour 
0|, which disperse with a particular nesting wave- vector 
corresponding to the allowed spin states of the warped 
contour In contrast, the cross-sectional contours of 
the SSB in Bi2Se3 are invariably convex, and together 
with weak scattering amplitude for non time-reversed 
states could explain the lack of QPI, despite the observa- 
tion of local impurities and defects in this system [l6|, [ItJ . 
However, these defects do affect the local surface poten- 
tial and therefore should also affect the Local Density of 
Sates (LDOS), exhibiting resonance phenomena similar 
to d-wave superconductors [10] and graphene 11]. 



In this paper we present a detailed study of defects- 
induced resonance states on the surface of the 3D topo- 
logical insulator Bi2Se3. Different defects are studied, 
all showing evidence for the creation of low-energy res- 
onance states which vary in strength but do not locally 
destroy the Dirac point. We further show that the main 
features of the data are captured by a simple extension 
of theoretical analysis of point-impurity scattering in the 
presence of Dirac dispersion [1], thus stressing the ro- 
bustness of the SSB. In particular we present high reso- 
lution topography and spectroscopy studies of resonance 
states associated with the frequently observed triangular 
defects, while also making a direct connection between 
their occurrence and the samples bulk doping. 

Single crystals of Bi2Se3 were grown by slow cooling a 
binary melt as described elsewhere [l4| , yielding samples 
with typical dimensions of 2 x 1 x 0.1 mm 3 . The trigonal 
c-axis is perpendicular to the cleavage plane of the crys- 
tals. Low magnetic field Hall measurements of similar 
undoped samples indicated n-type doping of ~ 1 x 10 19 
carriers/cm 3 which is a result of Se deficiency [ill, fl4j . 
In this paper we also show results on Sb-doped sam- 
ples. These lower carrier density samples (about 100 
times lower), with much reduced triangular defects, were 
obtained by slow cooling a ternary melt containing pro- 
gressively more Sb. Although Sb is isovalent with Bi, its 
substitution acts to control the defect density in the bulk 
crystals, reducing the bulk carrier density [15 ]. 

Samples were cleaved in vacuum of better than 5 x 
10 -10 Torr, and quickly lowered to the ~ 8 K section 
of the microscope, where cryopumping ensures that the 
surface remains free from adsorbates for weeks or longer. 
Topography scans were taken at several bias voltages 
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FIG. 1: 1600xl000A 2 topograph of the surface of undoped 
Bi2Se3- Most impurities are triangular (typical one marked), 
pointing in the same direction, and are a consequence of a sub- 
surface Se vacancies [16|, LL7j] . Inset shows a three dimensional 
rendering of a triangular defect. 



and setpoint currents (typically -f lOOmV and 60pA). Fig- 
ure Q] shows a 1600 x 1000 A 2 topograph of the unoccupied 
states at the surface of cleaved undoped Bi2Se3 crystal. 
Similar to previous studies of undoped Bi2Se3 crystals 
[ill [l7| , we observe a large concentration of triangular 
defects, uniquely oriented, with similar size but different 
brightness. High resolution topographs of the brightest 
triangular defects (see inset of Fig. [1]) revealed a structure 
of six "bumps," each about ~ 0.1 A high, separated by a 
distance of 8.2±0.2A, corresponding to twice the hexag- 
onal lattice constant of Bi2Se3. Less intense triangles 
have the same spatial structure but weaker in strength. 
Assuming Se vacancies as the main source of defects in 
our samples [14[, we can interpret the triangular struc- 
tures as a response to the Se vacancies in different layers 
beneath the surface. As we will argue below, this in- 
terpretation is further supported by an earlier study of 
Urazhdin et at. [18|, where excess Bi was introduced to 
Bi2Se3 crystals, producing clover-shaped defects in the 
occupied-states topographs, with the same dimensions 
as the triangular defects. Detailed calculations of these 
defects [18| suggested that they result from antisite dop- 
ing where a Bi atom replaces a Se atom in the bottom 
Se layer of the top quintuple stack. The fact that reso- 
nance states can be identified in that paper (see Fig. 3a in 
[H|), with opposite charge compared to our data below, 
further supports our proposal. Counting ~ 60 triangu- 
lar defects in this topograph, and assuming that they all 
come from within the first quintuple layer, each donating 
two electrons per defects, we find a bulk carrier density of 
~ 6 x 10 18 carriers/cm 3 similar to the one deduced from 
Hall effect. The much reduced triangular defects found 
for the Sb-doped samples, for which also the bulk carrier 
density is low, supports this scenario. 

Figure [2]depicts a topograph of a region containing sev- 
eral triangular defects, together with LDOS at three ener- 



FIG. 2: A topograph (a) , and selected energy maps em- 
phasizing two triangular impurities with different intensity, 
while both are the same size, bright triangle corresponds to 
a stronger protrusion than the faint one. LDOS at 90 meV 
(b), —180 meV (c) (points 1, 2, and 3 mark the center of de- 
fect, side-resonance, and background respectively) , and —400 
meV (d) are also shown. 



gies: close to the bulk conduction band (BCB) (Fig. [2)), 
in the middle of the exposed Dirac SSB (Fig. [2t), and 
below the Dirac point close to the bulk valence band 
(BVB) (Fig. [2]i). In analyzing these maps (and many 
similar maps we studied), we concentrate on two trian- 
gular defects with different brightness, corresponding to 
a different amount of protrusion. The first thing to notice 
that while LDOS in the presence of the BCB is clearly en- 
hanced on the defect site (Fig. (2)3), the LDOS in the pres- 
ence of the BVB is suppressed (Fig. [2]i). This indicates 
that the defect site is positively charged, which again is 
confirmed by the fact that our undoped Bi2Se3 is n-type. 
However, the most interesting behavior is found within 
the exposed Dirac-SSB. As is seen in Fig. [2^, three clear 
LDOS peaks around the triangular defects are formed, 
whose intensity strongly depend on energy. Although 
the lack of QPI implies a suppression of backscattering, 
the residual scattering of other momentum-transfer wave- 
vectors add up to a give rise to a resonance state [1]. 

The Sb-doped Bi2Se3 samples show reduced density 
of triangular defects, but exhibit other types of defects. 
Fig. [3] shows a topograph of a typical defect (presumably 
Sb-cluster) and LDOS maps of the same region. Here 
again, for energies overlapping with the BCB and BVB 
the LDOS at the defect site and around it is depressed, 
while a strong peak, now spread at and around the impu- 
rity, appear in the exposed region of Dirac-SSB (Fig.[3t). 

Figure 2^ shows LDOS on the triangular impurities 
("bright" and "faint"), with a LDOS curve away from 
impurities shown for reference (2p| . Fig. [4)3 shows the dif- 
ference in LDOS between those curves, showing a broad 
resonance peak that is located at positive bias as com- 
pared to the Dirac point. Both resonances are rather sim- 
ilar in shape, and we emphasize here that similar shape 
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FIG. 3: A topograph (a) , and selected energy maps empha- 
sizing an extended impurity in Sb-doped Bi2Se3. LDOS at 75 
meV (b), —225 meV (c) (relevant to Fig. El points 1, 2, and 
3 mark the center of the defect, the edge of its extent, and 
background respectively) , and —400 meV (d) are shown. 



and width have been measured for many triangular de- 
fects with a variety of "brightness" on several different 
samples. Fig. 0t shows LDOS on top of the defect shown 
in Fig. [3] on a Sb-doped crystal, together with the LDOS 
away from the impurity. Fig. 3]i shows the difference 
curve for these two spectra, resulting in a much sharper 
and more pronounced resonance. 
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FIG. 4: Emergence of a resonance near the Dirac point for 
two different impurities; a) Triangular defects correspond- 
ing to the "bright" and "faint" of Fig. [2) with the difference 
LDOS for both defects shown in (b); c) Extended defect in 
Sb-doped Bi2Se3 (Fig. [3}, and the difference LDOS shown in 
(d). Dashed lines mark the location of the Dirac point. 

The general problem of a impurity-potential scattering 
on the surface of a TI has been recently studied by Biswas 
and Balatsky Considering a two-dimensional Dirac 
gas interacting with a (^-function-potential, they found 
that a resonance state is formed around any localized de- 
fect; the Dirac point is not destroyed by the scattering; 
and LDOS around the impurity falls of as 1/r 2 . Further- 



more, the shape of the resonance is not Lorentzian and 
its width is determined by the strength of the potential 
which also sets the position of the resonance on the en- 
ergy axis. The stronger the potential is the closer the 
peak gets to the Dirac point and the narrower it gets. 
While neither the triangular defects nor the Sb-cluster, 
can be approximated by a ^-function potential, we as- 
sume that the net effect on the SSB may still be similar, 
especially if we look at the LDOS slightly away from the 
center of the defect. However, to account for the ex- 
tended nature of the defect we modify the calculations 
of ref. but with a step-potential replacing the local 
^-function potential. The step-potential is parametrized 
with a height V — Vb, extending over a radius R < Rq, 
and zero everywhere else. Using a step-potential, the 
strength of the impurity is therefore characterized by 
the product VqRq, rather than the potential height it- 
self. The extent of the potential (Rq) is further checked 
against the topography and the LDOS in the exposed 
part of the Dirac band. In particular, since the center 
of the triangular impurity shows suppressed LDOS, the 
"on-site" LDOS for these impurities are calculated as an 
average between the LDOS at the center and the LDOS 
at the three peaks surrounding it (shown in Fig. [5^). 
Fig. [SJd-c shows typical fits to the measured LDOS dif- 
ference for both, triangular and Sb-cluster defects. 

The first thing that we notice is that the quality of 
the fit is visibly better for a triangular defect (Fig. [SJd 
compared to the extended defect on the surface of an Sb- 
doped sample (Fig. [5]i). We believe that this is due to 
the fact that the triangular defects are more compact in 
space and are therefore better described as a local pertur- 
bation. The Sb-cluster on the other hand is a significantly 
extended object whose effect might not be described as 
a simple potential interaction with SSB. For example, 
in a recent study of an extended step on the surface of 
Bi2Te3 we have shown that a "bound- state" is formed 
as a consequence of interference in the bulk due to the 
finite penetration depth of the SSB [21[. Specifically, an 
important feature of the theory is the protection of the 
Dirac point against the perturbation due to the impurity. 
Here, as the bare LDOS (far away from any impurity) is 
subtracted from the LDOS that is affected by the impu- 
rity, and assuming that the background incoherent LDOS 
is uniform in space and only weakly energy dependent, 
the zero LDOS should occur at the Dirac point. While 
this feature is satisfied for the triangle impurities, the Sb- 
cluster impurity shows excess scattering near the Dirac 
point, shifting the zero of the difference-LDOS (Sp) by 
about 30 meV. However, as is seen in Fig. [SJi away from 
the center of the defect the fit becomes much better, in- 
cluding the location of the Dirac point (zero crossing). 
Therefore, the overall fit of our model to the extended 
defect data captures the shape of the resonance, and the 
fact that being a strong impurity it is sharper and peaks 
closer to the Dirac point. 
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Focusing on the triangular impurities, we already men- 
tioned that they appear with a variety of strengths re- 
lated to their different protrusion heights (see Fig. [TJ) A 
surprising result, which is revealed when we examine the 
LDOS difference of the "bright" and "faint" triangular 
defects, is that they have the same scattering strength 
as both positions of the associated resonances and their 
peaks widths are the same within the error-bars. This is 
in a good agreement with the suggestion that all triangle 
impurities have similar origin. 
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FIG. 5: Theoretical fits to resonance amplitudes, a) raw 
LDOS difference data for a "bright" triangle where (1) rep- 
resents the middle of the triangle and (2) represents the side 
peaks (see Fig. [2k). The averaged data is marked with Av.; 
b) Averaged data from (a), fitted to theory (see text) with 
Vb = 3.17 eV, R = 15 A; c) Fit to extended defect on an Sb- 
doped sample, at center of defect (point (1) in Fig. [3]) with 
Vb = 0.71 eV, R = 64A; d) Same parameters as for (c), but 
away from center of defect (point (2) in Fig. [3} (see text). 
Dashed lines mark the location of the Dirac point. 



FIG. 6: LDOS as a function of distance from the center of the 
impurity. Dashed line is a fit to 1/r for the dark points (after 
~ 65 A) and dash-dotted line is a fit to 1/r 2 . Inset shows the 
contours of decay of the LDOS around the impurity center 
(•) at 65A, lOOA, and 150A. 



rise to resonant states. Comparing the various defect res- 
onances to each other and with theory, we conclude that 
a simple scattering mechanism in the presence of a Dirac- 
type 2D quasiparticles give a full account to the observed 
data. While theoretical curves fitted the data very well 
for the case of the triangular defects, we find that for 
stronger, more extended defects, some modifications are 
needed that may include the finite penetration depth of 
the surface state. 
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06NA25396 and UCOP TR-027. 



Finally we show in Fig. [6] the dependence of the LDOS 
at the energy of the resonance of the broad scatterer (Sb- 
cluster) on the distance away from the impurity center. 
Fitting the decay as a function of distance, we find it 
hard to determine the exact power law decay. In Fig. [6] 
we show an averaged decay for the Sb-cluster defect, fit- 
ted with both, 1/r and 1/r 2 . While the 1/r decay fits the 
data better all the way to the extent of the defect, the two 
power laws, and in fact any power law in between may 
still be a good description of the data, especially consid- 
ering the large characteristic length for Bi2Se3 surface 
state (A ~ 30 A |9[). In particular the data rules out an 
exponential decay of the LDOS resonance with A, which 
is another confirmation of the two-dimensional nature of 
the band giving rise to the resonance state. 

In conclusion, we identified the resonance states associ- 
ated with the triangular impurities in Bi2Se3. We found 
that independent of their observed intensity, these im- 
purities induce resonances of similar shape and strength, 
thus can be considered as a characteristic of Bi2Se3 (and 
possibly also L^Tes) spectroscopy. We further showed 
that other type of defects, extended in nature, also give 
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